NACA TN 2666 9 £0 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 






TECHNICAL NOTE 2666 


TWO-DIMENSIONAL SUBSONIC FLOW PAST ELLIPTIC 
CYLINDER BY THE VARIATIONAL METHOD 
By G. V. R. Rao 
New York University 



Washington 
March 1952 



TECH LIBRARY KAFR, MM 



TECH LIBRARY KAFB, NM 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 




TECHNICAL NOTE 2 666 


TWO-DIMENSIONAL SUBSONIC FLOW PAST ELLIPTIC 
CYLINDER BY THE VARIATIONAL METHOD 
By G. V. R. Rao 


SUMMARY 


A method of solution is presented for compressible fluid flow past 
an elliptic cylinder by means of the variational method. The solution 
is obtained- as a function of thickness ratio and free -stream Mach 
number. Numerical examples are carried out for several thickness ratios 
and Mach numbers and the results are compared with those obtained by 
other methods. It is seen that the variational method yields good 
results for flow past a thick body at a low Mach number as well as for 
flow past a thin body at a high Mach number. 


INTRODUCTION 


In recent years the progress made in theory and practice of high- 
speed aerodynamics has made all the more important the fact that one 
should try to solve the equations of motion more accurately. Owing to 
the difficulty of the nonlinearity of the equations, this problem has 
been attacked always by approximations. 

For slender bodies there is the solution given in terms of the 
thickness parameter, usually referred to as the Prandtl-Glauert method. 
On the other hand the method of Rayleigh-Janzen gives the solution in 
powers of Mach number and can be used in the case of thick bodies. This 
method was applied by Hooker (reference 1) to the study of compressible 
fluid flow past elliptic cylinders. But owing to the necessity of 
expanding a certain function in the analysis. Hooker's method cannot 
suitably be applied for small fineness ratios. Kaplan by using Poggi's 
method has investigated this problem (reference 2). Although his result 
is expressed in f init e terms, some infinite series must be used in the 
course of the method. Moreover, Kaplan's solution deals with the 
velocity and pressure distributions on the surface of the cylinder only 
and not in the interior of the domain. Imai and Aihara (reference 3 ) 
also used the Rayleigh-Janzen method and obtained the flow past an 
elliptic cylinder in a more general manner. The use of the Rayleigh- 
Janzen method would involve formidable computation beyond the second 
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approximation. For tills reason this method is not very suitable . f or 
flows at high Mach numbers. Perl (reference 4) used the method of 
expressing the equations of motion in terms of the streamline curvature 
and then integrating them. 

In this paper is presented a solution for compressible fluid flow 
past an el l iptic cylinder by the direct method following the calculus 
of variations. The formulation of the variational principle in hydro- 
dynamics was first due to Hargreaves (reference 5) ‘who has shown that 
the integral to be used in the variational, principle is a linear function 
of pressure . For steady irrotational flows, Bateman (reference 6) has 
shown that the pressure energy is an extremum. He gives the Integrand 
as a function of the density and the velocity potential. In an earlier 
paper (reference 7) Bateman gives the integrand, as a function of the 
velocity and local velocity of sound. Bateman’s original integral, 
however, can not be directly used for fluid which occupies an infinite 
domain. Braun (reference 8) used the same integral without modification 
and attempted to obtain an approximate solution for uniform steady flow 
past a circular cylinder by the use of Rayleigh-Ritz method and his 
method leads to erroneous results. Wang (reference 9) modified the 
integral for such cases and it is fundamentally this method that was 
used by the author in this report for the study of flow past elliptic 
cylinders . 

An elliptic coordinate system is used and in the set of simultaneous 
equations given in the appendix the coefficients of the terms contain 
as parameters the thickness ratio of the ellipse and the free-Btream 
Mach number. Thus a method of solution for the flow past an elliptic 
cylinder of any fineness ratio is indicated. The results are compared 
with those obtained by Hooker, Kaplan, and Perl. 

The present work, carried out at the Daniel Guggenheim School of 
Aeronautics, College of Engineering, Hew York University, is a major 
part of a thesis submitted to the Graduate Division for the partial 
fulfillment of the degree of Doctor of Engineering Science and has been 
made available to the NACA for publication because of its general 
interest. 

The author wishes gratefully to acknowledge the valuable sugges- 
tions and guidance given to h im by Dr. C. T. Wang in the preparation of 
the thesis. The author also wishes to express his appreciation to 
Professor F. K. Teichmann for the facilities afforded to him. 



NASA TN 2666 


3 


LIST OF SYMBOIS 


a 

C P 

g = q m 2 - q 2 

I 

M 

P 

q 

U 

u 

v 

x ^y 


velocity of sound 
pressure coefficient 

variational integral 
Mach number ( q/a) 
pressure 
local velocity 
maxi-mum velocity 

velocity of uniform parallel flow at infinity 
velocity component in x-direction 
velocity component in y-direction 
rectangular coordinates 


z = x + iy = c cosh (I + in) 

7 ratio of specific heats, 1.405 for air 

T] eccentric angle 

• elliptic coordinates 
p density 

0 velocity potential 

80 admissible variation of 0 

Subscripts: 

1 incompressible flow 
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o conditions^ at infinity 

x,y,^,T| respective partial derivatives 


EQUATIONS OF MOTION FOR C02-URESSIBLE FLUID FLOW 


The various equations to be satisfied for steady, irrotational flow 
of a compressible fluid are given below. The continuity equation in 
two dimensions is 


J^-(pu) + ^-(pv) = 0 
ox dy 


( 1 ) 


For the flow to be irrotational one must satisfy another relation 


5v chi _ q 
dx dy 


( 2 ) 


For isentropic changes of state, the equation of state is 


p = Constant X 


The Euler equations of motion in the integrated form will then be 


a 


2 




(3) 


where 7 is the ratio of specific heats and a is the velocity of 
sound. The subscript o refers to the conditions at infinity. Since 
in the variational principle, as explained later, a value of 7=2 
will be used, the above equation of state is to be replaced by 


p = A + Bp^ 

where A and B are constants. The significance of such an equation 
and the justification of using 7=2 are given in reference 9* This 
equation of state in the differential form is 
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— = B7P^ - ^ = 


dp 


= a c 


( 4 ) 


One substitutes equation (4) into equation ( 3 ) and introduces 


u 



v 


50 

dy 


Thus q 2 = 0 X 2 + 0 y 2 . By this introduction of 0, the irrotationality 

condition (equation ( 3 )) is satisfied a priori. Equation ( 3 ) will now 
take the form 


B7P 7 " 1 + + (/) 



Since the right-hand side of this relation is a constant, it can he 

replaced hy 21_Z c^ 2 where q^ 2 is dependent only upon the conditions 

at infinity . Then 


B7P’ - - 1 = 2L §^(<J m 2 - fi x 2 - fiy 2 ) (5) 



Substituting this value of p into the equation of cont inuit y (l), one 
obtains the differential equation for the velocity potential 
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The "boundary conditions on the velocity potential are as follows: 
When the flow is subsonic, the disturbances caused by the body would 
vanish at infinity. Thus the velocity at infinity would be that of the 
uniform parallel flow. The boundary of the body would be a streamline 
and so the normal derivative of 0 is zero on the body. Expressed 
mathematical 1 y, the boundary conditions are 




M = u 

dx 


at infinity 

► 


M = 0 
dn 


on the body J 


(7) 


VARIATIONAL IRINCIPLE 


The formulation of a variational principle for the motion of a 
nonviscous compressible fluid is given by Bateman (reference 7)* 

He indicates the associate'd integral in the calculus of variations is 
of the type 


I = 


G(q) dx dy 


Further, he suggests that for adiabatic conditions the appropriate form 
of G is 


G(q) 




( 8 ) 


2 7-1 2 

where qQ is the velocity at infinity. Replacing a Q + ■ — — q G 

by — ~ — q^ 2 results in 


G(q) 


= (7 - 1 ) 7 “! (c 


% 



7 

7-1 
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Since the constant coefficient is immaterial in a variational principle, 
Bateman’s integral can he written as 

1 = J (\ 2 ~ q 2 ) 7 ' 1 dx dy (9) 

A proof that the integral, equation (9)j gives equation (6) as its 
Euler equation can he found in reference 9* 

For irrotational flow q 2 may he replaced hy 0 X 2 + 0 i then 

w 

I = yY (q m 2 - 0 X 2 - 0y 2 ) 7-1 dx dy (9a) 

Now 0 is varied so that the variation 50 is zero on the houndary 
which is not a streamline^ 81 = 0 becomes 

1 

(Sn 2 “ ^x 2 “ 0y 2 ) 7-1 (0x 6 0x + 0y 5 0y) dx dy = 0 (10) 

Equation (10) can he modified as follows: 




dx dy = 0 


( 11 ) 
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With the aid of the divergence theorem, the first integral can he trans- 
formed into a line integral taken on the boundary. This is 



The line integral taken on the body contour will be zero since ^ = 0 

8n 

from boundary condition (7)- The variation 50 must be so chosen as 
to approach zero on the other boundary. For flow in an infinite domain 
a correction must be made so as to insure the vanishing of the line 
integral as was carried out in reference 9* Although 50 approaches 
zero at infinity the line integral taken on the infinite boundary may 
be of finite value. This should be subtracted from integral (9)> so 
that from equation (ll) the following expression remains 



(13) 


Since 50 is any admissible variation, it follows that to satisfy this 
relation equation (6) must be satisfied in the whole domain. 


ELLIPTIC COORDINATES 


For solving two-dimensional flow past elliptic cylinders, it is 
convenient to use elliptic coordinates and transform the integral, 
equation (9) , into the elliptic coordinate system. 

Take | and q defined by the equation 

z = x + iy = c cosh (| + irj) 
where c is a read, constant. 
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The above relation after separation of real and imaginary parts 
would give 


X = C COSb | COS T) 

y = c sinh £ sin t] 


One takes all values from zero to infinity for I and values from zero 
to 2rt for T). The curves I = Constant and t) = Constant will be 
the directions of the coordinate axes and are confocal ellipses and 
hyperbolas^ respectively. These are shown in figure 1'. The angle T] 
is usually called the "eccentric angle." 

If a and b are the semimajor and semi-minor axes of the ellipse 
given by £ = § , then 


a = c cosh £ 
s o 

b = c sinh £ 

o 

_ p p 0 

It follows that a - b = c , and 


a + b 
a - b 



The ratio b/a 

the function e 
computation. 


is the thickness ratio of the ellipse and for simplicity 




o _ 


a + b 


is termed t in the author's discussion and 


The modulus of transformation is given by 


d( I + it)) 
dz 


1 

c(cosh 2 i - cos 2 ^) 1 / 2 


m 
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If f is tbs velocity potential in elliptic coordinates, the velocity 
components normal and tangential to an ellipse are 

h K 

and - — ; ; — — , respectively. Hence 


:^cosh 2 | - cos 2 t 


:^cosh^| - cos 


q 2 = 


<t>i + K 


c^cosh^I - cos 2 t]) 


(15) 


Also the elemental area in the elliptic coordinate system is 


da = 


d£ dq 

c 2 (cosh^| - cos^q) 


(16) 


APPLICATION OF RAYLEIGH-RITZ METHOD 


After transforming the b oundary -value problem into a variational 
problem, the Rayleigh-Ritz method can be used to obtain close approxima- 
tions to the solution. Using relations obtained in the previous section, 
variational integral (9) in the elliptic coordinate system would be, 
with a value of 2 for 7, 


c^cosh^i - cos 

and it is proved in reference 6 that for subsonic flow the integral I 
is a maximum for all admissible variations of 0. 

The Rayleigh-Ritz method is a direct method and consists of 
approximating p by a finite sequence of functions. For 0 in the 
integral (equation (17) )> substitute 

<f> = H a k 0, (18) 

k=l 




c 2 (cosh 2 | - cos 2 q) d£ dq 


(17) 



\ 



NACA TN 2666 


11 


with. 0^ satisfying the boundary conditions and with as undeter- 

mined parameters. 

Since the required solution for 0 makes the integral a maximum, 
one may transform the variational problem into an ordinary maximizing 
problem. Thus one can obtain a system of simultaneous equations by 
taking partial derivatives with respect to the various constants: 





= 0 


and so forth. Solving these simultaneous equations, the constants can 
be obtained. Then 


0 = a 1 0 1 + a*p 02 + . . . (l8a) 

makes the integral (equation ( 17 )) stationary and is a solution of the 
problem. 

In the case where more than one set of values for these constants 
is obtained the solution which makes the integral a maximum is to be 
taken as the required solution. 


EXTENSION TO HEINITE DOMAIN 


In applying the variational integral, equation { 9 ), to the present 
case when the domain extends to infinity, the value of the integral 
becomes infinite. But this does not prevent one from obtaining a solu- 
tion by maximizing or minimizing such an integral. The difficulty 
would arise only in the case when this integral has to be evaluated. 

There is another feature in connection with line integral (12) 
which should be scrutinized carefully since it is taken around an 
infinite boundary. This should be zero to satisfy the conditions upon 
which the variational problem is formulated. Even though 50 approaches 
zero at infinity, the product 50 ds may contain certain terms of zero 
power in r. Then line integral (12) will have a finite value and will 
have to be accounted for by subtracting it from the integral of 
equation ( 9 ) • 
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Consider line integral (12) around an ellipse given by a large 
value of | and then let | approach infinity. One obtains for 





2 

- COS T) 


j 


with 0 assumed as given in a later section, equation (23) > at large 
distances 


where 


°f— 2 

\cosh i) 


q 2 a u 2 + ; 

\cosh^| / 

means terms approaching zero to the order of l/cosh% 


and higher as | — 

In the line integral, ^ is in the direction of the inward normal. 

on 


Therefore 


^ = h 

Sn SI 


and the elemental length 


ds = g- dT) 


where h is the modulus of transformation. Hence 


Mi, Jj, 

Sn S| 

and with 0 as in equation (23), at large values of |, 


Sn 




ds = U(a + b) sinh | cos T) + 0( 


o 

cosh I , 


dt) 


(19) 
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In equation (23), the assumption is that the flow at infinity is pre- 
scribed and is equal to the incompressible flow. Hence, the -variation 
60 for large values of I will be the same as the difference between 
the assumed 0 and the incompressible velocity potential. Thus 


80 = Ap! 


£22_JL + A £SS^L + of— i— l 

cosh I ^ cosh I \cosh3 ij 


( 20 ) 


Substituting these into line integral (12) 



After integrating and letting 5 — >°°, the above expression becomes 


- ^(q m 2 - U 2 )u(a + b) ^A n + l «A 13 ^ (21) 

This is the term which should be subtracted from the integral of 
equation ( 17 ) • 


FLOW PAST AN ELLIPTIC CYLINDER 


The velocity potential for incompressible flow past an elliptic 
cylinder defined by I = i Q is 


0^ = U(a + b) cosh (I - £ 0 ) cos q 


( 22 ) 
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where a and b are semiaxes of "the ellipse. Tbis can be readily 
decked to satisfy Laplace’s equation in elliptic coordinates 



To investigate the two-dimensional flow of a compressible fluid 
past an elliptic cylinder given by £ = Iq, one may assume 0 in the 
following form: 


a a » cos 

0 = 01 + A U coshTT 


So' 


"^33 


3 

COS J T] 


cosh^d - 5 n ) 


(23) 


where 0 ^ is the potential for incompressible flow and .A^-^ A-j^j ^ 31 * 
and A 23 are undetermined parameters. 

In the above assumed potential (equation (23)) all the additional 
terms approach zero at infinity to the order of l/cosh |. This leaves 
only the potential of the incompressible flow which is the boundary 
condition at infinity. Also it can be seen that for the additional 
terms the derivatives with respect to i vanish at | = £ Q . Thus the 
normal velocity component on the body due to the additional terms is 
zero. The assumed variation (equation (23)) therefore satisfies the 
required boundary conditions. If the flow were not symmetrical, one 
would have to add to expression ( 23 ) terms contai nin g sin 3 , sin3ri, 
and so forth. 


The assumed 0 function (equation (23)) is now substituted into 
the variational integral, equation (9)» To cover the whole domain 
exterior to the ellipse given by | = £ 0 , the limits for integration 
are from 0 to 2it for q and from to co for | . By differ- 

entiating under the. integral sign, one obtains the following four 
s imult aneous equations: 


dl 

SaH 


0 



dl 


3a 


= 0 


13 





> 


j 


(2b) 
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Integrating these equations, four algebraic simultaneous equations in 
the parameters Ay), A^ 3 , A^, and A33 are finally obtained and are 

given in appendix A. The equations given in appendix A are divided 
into two parts. The linear terms in the parameters are on the left- 
hand side of the equations. All the non lin ear terns are grouped 
together and called K-j_, Kg, K3, and Kjy These are given at the 
right-hand side of the equations, together with the constant terms. 


NUMERICAL EXAMPLE 


The flow past an elliptic cylinder of thickness ratio b/a = 0*5 
when the Mach number of uniform flow at infinity is 0.5 is investigated 
and the method used in solving the equations is illustrated here. 

The value of the major .axis of the ellipse is taken as a unit of 
dimension, with b/a = 0.5* Then 


t = = 3 

a - b 


Divide throughout the equations of appendix A by the velocity of sound 
at infinity, so as to make the coefficients nondimenBional. Substituting 

these values of u/a 0 , q m 2 , and t in the four simultaneous equations 
of appendix A and evaluating all the coefficients, one obtains the 
following system of equations: 

1.99301954^* + 1.67206143^3 ' + i.62409948a 31 * + 1.32770102A 33 * = 
0.09169413 + 2^5 V 


I.67206426A3J/ + 2.28737268A33 * + i. 3832566 a 3 3_* + l.7l470758A 33 * = 
0.03621830 + k 2 * 

1.62409948A3JL* + 1.38325662A 13 * + i.96462204a 3] _ * + l.59203403A 33 * = 
0.07905474 + K 3 * 
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1 • 32770102A 11 * + 1.71470758a 13 * + 1.59203403A 31 ' + 1 . 78968390 A 33 ’ 


0.03812186 



where K^*, Kg’ , K 3 ’, and K ^ 1 are the proper values of K^ Kg, 

K 3 , and K 4 of appendix A with the values of t, u/a 0 , and q m ^ 
substituted and Ain*, A 13 1 , and so forth in the above set of equations 
are the assumed coefficients of appendix A, each divided by the velocity 
of sound at infinity. 

Wow one may make use of Crout's method (reference 10) and solve 
the above set of equations in terms of the right -hand- side members of 
each equation. Since the above matrix is symmetrical, it is quite easy 
to obtain the following expressions: 


a^jl* = 4.670656231! - 3. 95643729kg - 4. 79816972X3 + 4. 5396629X4 
Ai 3 * = -3-95643629X! + 4.99700942kg + 4.49078519X3 - 5.84736055X4 

a 31 ' = -4.79816982k! + 4.49078533X2 + 6.86307052X3 - 6.814818730X4 

A 33 ' = 4.59396635X-L - 5.84736o64kg - 6.814818723X3 + 8.84494015X4 

where 

X-i = 0.09169413 + -i— K n * 

1 2.25 1 

Xg = 0.03621830 + Kg’ 

X 3 = 0.07905474 + k 3 « 

X 4 = 0.03812186 + K 4 ’ 
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For the computation of A-y ' , first one puts all the nonl in ear 
terms as zero and obtains the values of the coefficients. Next these 
are substituted in the complete equations and the values of the coef- 
ficients are again solved for. This procedure is repeated and the above 
set of equations is solved by successive iterations, until the coef- 
ficients iterate themselves within a desired accuracy. The following 
are the values obtained: 


A-jjl* = 0.08335840 

a 13 * = - 0.05016832 
A 31 * = 0.00714746 
A^’ = 0.00157309 


(25) 


These values are substituted in the expression for the velocity potential, 
equation (23). The velocity at a point defined by the eccentric angle T) 
on the elliptic cylinder is given by 


q = 


C^COSh % 0 - COB^tl J - 1 -/ 2 ^ 


The value of is obtained from 


t - e 2? ° - 3 

With the value of 0 substituted in the above equation and after 
dividing through by the velocity U at infinity, q/u becomes 


1 = 
U 


c^cosh 2 | 0 - cos 


W 1 


(»•?) 


An’ 

sin Ti + — - — sin ti + 
Mo 


3A-,.’ 


131 2 _ „ . ^ 31 _* 


Mo 


cos T) sin q + 


Mo 


3^-33’ 2 

■ cos t] sin tj 

% 


sin tj + 


(26) 
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With the obtained, values (equations (25)) of the coefficients, b/a = 0.5, 
and Mq = 0.5, the above equation would be 


U = 1,1 ^ 701 P^ (i.68101172 - 0.2915H38 cob 2 ti) sin t \ (27) 

(3 ' cos2 T 


These values of q/U are computed for various angles n and are 
given in table 1. The results are also shown plotted in figure 2. 


From the pres sure -density relation for adiabatic flow one can 
obtain (reference 11) the following relation between the local pressure 
and pressure at infinity: 


P - Po 

K ° 2 



(28) 


Note that 7 = 2 was used in the process of obtaining a solution for 0. 
It is possible to use 7 = 1.5 or 7 = 1*33> but it is seen that the 
cha n ge in 7 does not alter the computed values of the velocities 
appreciably. 


Even though 7=2 is used for getting the velocities it is 
proposed to use 7 = l.b05 in computing pressures. Then 


C„ = 


P - Po 


P !L tt 2 


0.7025V 


1 + 0.2025V 


[‘ - ( 8)1 


^> 1.405 

oT5o5 


(2 9) 


The pressure computed as' above is shown in figure 3, as a pressure 
coefficient, against T] and is compared with the pressure distributions 
obtained by the Prandtl-Glauert method and the Karman-Tsien method 
(reference 12). 

By substituting t = ll/9 in equation (2b), the flow past an 
elliptic cylinder of thickness ratio 0.1 is computed when Mq = 0.8, 
and the results are shown in figures b and 5> 
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Using Poggi’s method Kaplan (reference 2) has calculated the flow 
past elliptic cylinders of thickness ratios 0.5 and 0.1. These results 
are shown in tables 1 and 2 and are plotted in figures 2 and 4 for 
comparison. 

By letting t approach infinity the flow past a circular cylinder 
is obtained and the result is shown in table 3 and figure 6. Imai 
(reference 13) has calculated the subsonic flow past a circular cylinder 
to the sixth power of the stream Mach number. These results are shown 
in table 3 and figure 6. The results obtained by Wang (reference 9) 
show a higher velocity throughout. This is due to the fact that Wang 
used six terms in the series assumed for the velocity potential. In the 
present investigation only four terms of the series are -used and hence 
lower values for the velocities are obtained. 


CONCLUDING DISCUSSION 


The velocities obtained by the variational method show more devia- 
tion from the incompressible flow than indicated by Kaplan, Hooker, and 
other investigators. From the velocity distribution shown in figures 2 
and 4, it is seen that the variational method yields a. good result for 
flow past a thick body at a low Mach number, as well as for flow past a 
thin body at a high Mach number. There is close agreement with Perl’s 
results at both thickness ratios. 

The re stilts obtained by the variationhl method for the flow past an 
elliptic cylinder of 0.1 thickness ratio at a Mach number Mq = 0.8 are 
given in table 2. The maximum velocity occurs at the ends of the minor 
axis and Is given as q/u = 1.2255* With the aid of the well-known 
relation between velocity and the stream density for adiabatic flow, the 
maximum stream density is computed for the above velocity. A maximum 
stream density of 0.575 is obtained by such computation. Using the 
Prandtl-Busemann method, the second approx i mation carried out by 
Schmieden and Kawalki (reference 14) shows a maximum stream density of 
0.583- Further, Schmieden and Kawalki indicate that the value of stream 
density should natural 1 y coincide with 0.578 given by the sonic boundary. 
The present method gives results higher than those of Kaplan (refer- 
ence 2) and more in agreement with references 4 and 14. 

For the flow past a circular cylinder, the pressure coefficients 
obtained by the present method are between those of the Prandtl-Glauert 
and Karm^n-Tsien methods. 

For flow past elliptic cylinders the pressure coefficients obtained 
at T) = ji/ 2 are higher than those from the K^rmAn-Tsien method also. 
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From figures 3 and 5 it appears that the local Mach number has more 
effect than that indicated by the Karman-T s ie n method. In figure 7 the 
flow past a cylinder of 0.1 thickness ratio is shown at various Mach 
numbers. At higher Mach numbers the variational method gives results 
which differ considerably from those obtained by other methods. This 
deviation appears in Hooker’s results to some extent, as shown in 
figure 2. This may be due to the fact that the compressibility effect 
is far higher at points where the local velocity is nearer the critical 
value. The methods of Prandtl-Glauert and Kdrmdn-Tsien for obtaining 
the compressibility effect are independent of the critical Mach number 
of the flow. But in the variational, method the critical Mach number 
directly affects the coefficients in the velocity potential. In fact, 
at some value higher than the critical Mach number the coefficients do 
not converge on successive iterations. 

The point of zero pressure on the surface of a cylinder moves rear- 
ward when compressibility is taken into account. This is not the case 
when the Prandtl-Glauert and Karman-Tsien methods are used. The investi- 
gations of Kaplan and Hooker, however, show this deviation which is 
consistent with the experimental results. 


New York University 

New York, N. Y., August 1, 19^9 
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APPENDIX A 


EQUATIONS FOR FLOW PAST ELLIPTIC CYLINDER OF 

t - 1 


THICKNESS RATIO 


t + 


y ^ Mq = u/eq 


The equations for the flow past an elliptic cylinder of thickness 
ratio — y at Mq = u/a Q are as follows, where A, B, C, D, and 

v T J. 

so forth are certain defi ni te integrals contain in g the parameter t 
and are defined at the end of the appendix. 

First equation: 

A n|l.333333q m 2 - U 2 (6b + 6tK - 1. 18070817)] + 

AijjLov 2 - u 2 ^6b + 2.47817466 - 0.59035447 ^)] + 
A 31 [l.o6666667q m 2 - u 2 ( 72 tc + 24 tL - 1.54737540)] + 

|o.8q m 2 - U 2 ^72tC + 1.31446908 - 0.54340470 = 

3 -^(-f - 1 * - i-59m692 + | lose rrr) * (^r)% 


33 


~ t + 
Second equation: 


A u[q m 2 - U 2 ^6B + 2.47817466 - 0.59035447 + 

A 13 jl.33333333q m 2 - U^B + 2.85882167 - 0.02182589 - - 0.28250893 
a 31 [o-S'Jm 2 - + 1.10262516 - 0.57454014 ij] + 

A-33 [^n 2 " U 2 ^ 2 ^ + 1.^3683614 - 0.06345864 i - 0.27337794 ^] = 
3^- 1 + 1-78426410 - 0.19314718 |) + (~^) 2 K2 


3 2at / A 
t + 
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T hir d equation: 

Ajjt. 06666667^ - D 2 (72tC + 24tL - 1.54737540)] + 

*13 E • 8q m 2 - U 2 ^72tc + 1.10262516 - 0.57k5h0lk + 

jl • 21904693qjn 2 - U 2 (864D + 96tM - 0.29524866)j + 

A 33 jo.9l4285l9q m 2 - + 1.94284848 - 0.52783764 ij| = 

\ 2 

d 3 nif- 4 + 2051 - 1 - 2 'i'to 1 ‘39) + % 

Fourth equation: 

Au^.8q m 2 - U 2 (72tC + 1.31446908 - 0.54340470 ijj + 

A 13^m 2 " + 1-43683614 - 0.06345864 | - 0.27337794 + 

A 31 Jo.9l4385l9q m 2 - u 2 ^864d + 1.94284848 - 0-52783764 ijj .+ 

A 33 |l.02857l03q m 2 - u 2 ^864d + 2.06249040 - o.o846n88 J - 0.26469675 ^ 
u 3 t 2 +\ (-6B + o.2944i4o6 - 0.18223425 ^ % 
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The nonlinear terms K l> Kqj and K]^ are as follows: 


*1 " A ll 3 (3 a) + 3204 - 0.30476230) + A^^SD + 0.90624977 + 0.19262107 i - 0.12594008 A - 0.04971839 A) t 

t t^/ 

^3(2^6 0 + 2048tQ - 0.14097 1 108 + 0.09345024 i + 0.03760128 -1 + 0.01714176 -1) + 

A 33 3 (^^ 0 + 0.78763104 + 0.23067504 | - 0.08105616 A . 0.02979936 4jj + 

All 2 ® ~ glrtc + *41 . - 1.10000028) + A 11 2 a 13 ^96d + 1.77142776 - 0.442657 ^ + 

Ah £a 31 ^^ e + 384tn ' 0.46412160 + 0.05704704 | J + A-^A^ApS E + 1.68942240 - 0.35366256 ^ + 
*13^ ^T[(-24tC + 1.92993774 - 0.02555028 i - 0.27105478 ^ + 
a 13 2a 11^96d + 1.97916568 + 0.06166172 ^ - 0.25047667 Aj + 

A 13 2 a 3i(^^ e + s -°15190 + 0.17372192 | - 0.23433120 Aj + 

a 13^33 ^^ e + 2-72042070 + 0.69054222 i - 0.34588152 A - 0.14399928 A^ + 

a 31% E + 384tN - 1.16001280 - 0.17114112 + 

A 3i a Allf^ t F + !336tP - 0.14124920 + 0.12966912 | + 0.04409856 A^ + 
a 31 2a 13 (^|^ f * 3-95553600 - 0.20380672 i + 0.04409656 A^ + 
a 31 2a 33^^ 0 + 2-15970016 - 0.03670656 i + 0.11280384 A + 0.05142520 4jj + 
a 33 2 « E + 1-19026800 - 0.14974848 | - 0.24776064 + 

a 33 2a u(^^ f + 1-91459592 + 0.22063104 | - 0.18427392 + 

F + 2-71604124 + 0.77606208 i - 0.27906768 ^ - 0.13974606 Jj) + 

A 33 2A 3l(^5^ 0 + 2-25168768 + 0.42723072 | - 0.09931680 A + 0.05142528 ij) + 

A ll A 13 n t^l( J * 8tc + 3-54737568 - 1.03333332 + A^AjjU ^A<- 576 d + 192tM - 2.38094556) + 

AjjA^U t^T^' 5760 + 2-57143464 - 0-94285944 |) + AyA^o AjA(- 576 d + 2.57143464 - 0.94265944 |) + 

*13*33° t^r(-57fe + 2.87500632 - 0.01465512 i - 0.51822222 A^ + 

*31*33° E + 2.09204160 - 1.19231904 + AjiAjjA^^i B + 3.56932160 - 0.66923104 I) + 

*11*13*33 E + 3-79220400 + 0.27125568 | - 0.47769408 A^ + 

*11*31*33^5^ F + 3.61908864 - 0.47112192 | + 0.08819712 + 

*13*31*33 F + 3.89762544 + 0.46507776 | - 0.36557280 ^ 
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Kg = + 0.49047590 - 0.16538478 |j + 

a,, 3 ( 32D + 1.493750 + 0.53412758 — - 0.14506687 ri- - 0.14991404 A- - 0.03051809 -4-\ + 

\ t t 2 t 3 t*j 

a^i 3 a + 0.38134272 - 0.08972160 i + 0.03760128 i + 0.017110.76 -Aj + 

Aj3 3 ^^ G + 1.13567832 + 0.46610640 i - 0.11696238 ^ - 0.12707712 - 0.02943540 + 

A U 2 ' I t^\(' elttC + 1 - Tr 36878 k - 0.51666667 + A^^^d + 1.97916568 + 0.06166172 | - 0.250^7667 -ijj + 

a 11 2a 3i(^^ e + 1.781166080 - 0.331)61552 e + 3.896102 + 0.1356278*) | - 0.23881)706 —j + 

a 13% + 2.81)1)1)95 + 0.5223653!) | - 0.1)1226033 ^ - 0.15I160086 -Aj + 

a 13 2a ii^96d + 2.81)999937 + 0.62l)li9833 | - O.36896I106 i - 0.11)839629 + 

A 13 Sa 3 i(^^ E + 2-84542110 + 0.7li05!i070 i - 0.33572088 - 0.11)3083l»l) + 

A-^A^^ip E + 1). 1)9892738 + 1-75572360 i - 0.36951)288 -A - 0.1)3131528 -Aj + 

A 31 2n E + 1 -l 1 H2592o - 0.577U248 + A^l^ F + 1.955536 - 0.20380672 i + 0.01)1)09856 - 

* 2.27957512 + 0-347648 i - 0.1681)0701) ~j + 

A 31 2a 33 0 + 2.1)721)8768 + 0.50808576 | - 0.0836!)l60 ~ + 0.05142528 + 

a 33 2ij E + I-88874090 + 0.32336194 i - 0.37832 616 i - 0.14708412 + 

A 33 s A 11 ^2§p^ p + 2.71604124 + 0.77606208 i - 0.279067068 i - 0.13974606 -jj + 
a 33 2a 13(^T' F + 4.43507589 + 1-84144320 i - 0.28275813 ^ - 0.41674986 - 0.08774370 ^ + 

A 33 2a 3 i(^F^ g + 3.11637312 + 1.03692240 | - 0.16510608 i - 0.08204112 ij + 

A ll A 13 u + 3.85987548 - 0.05110056 i - 0-54210956 + 

AuA 31 n ^5^-576 d + 2.57143464 - 0.94285944 ±j + 

A U A 33 U t^l( _576l) + 2.87500632 - 0.01465512 | - 0.51822222 ~j + 
a 13 a 31 u t^r(- 576l> + 2.70833960 - 0.06203048 i - 0.52347038 + 

a 13 a 33 d :jrf^-S76D + *-37500434 + 0.97810890 i - 0.79421988 ~ - 0.30175410 ij + 

t • 

A 31 A 33 n ^l\"^ E + 2.38053600 - 0.29949 696 £ - 0.49552128 + 

A U A l3 A 3l(^t^ E + *-030300 + 0.3*744384 i - 0.46866240 + 

a 11 a 13 a 33 E + 5-44o84l4o + 1.38108444 i - 0.69176304 -A - 0.28799856 ij + 

a H a 31 a 33^3T^ f + 4.10292624 + 0.53652480 | - 0.35664768 ij + 

*13*31*33^^ F + 5-38417896 + 1.53128544 | - 0-56259792 ^ - 0.27990684 ij 


* 
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K3 = A^ 3 


l ll 3 (-jr e + 12845 


). 20232640 


.01901568 + 

l 13 3 E + 0.91722360 + 0. 


+ 0. 

8l92tn + 0. 


' t/ 

.23434728 A - 0.11444712 -A. - 0.04792344 

4 t 2 t 3/ 

H + 8l92tR + 0.18477056 + 


t 2 t , 

a 31 3 (^P ' h + 8l ^ 24E + 0.18477056 + 0.25878528 A + 0.08626176 -A + 0.04644864 -A + 0.0132 
A 33 3 r^ j p h + !• 03949568 4 0.43053768 i -'0:015120 A. + o.ooi8i44o i + 0.01327104 -^\ 4 
All 2 ° t^V" 288l) + 96tM * 1 - 1 3 <A 7 a 78) + A^A-^A^ E 4 1.78466080 - 0.33461552 i) + 
a 11 2a 31^^^ f + 15 36tP - 0.14124928 4 0. 12966912 A 4 0. 04409856 j-g'j + 
a U. Sa 33 F + 1 - 8 0954432 - 0.23556096 A 4 0.04409856 pj 4 


4 0.12966912 A 
0.23556096 | 4 0.04409856 

.35416980 - 0.03101524 A - 

n f xy t 

#***• 


3/66^ 

\ t* 

A ll 2A 3 lf^ F + 

, 2. /ij824 . 

a U^ 33 ^~ f + 1 

A-1 O^U _2S 4 _/. 2 88d 4 1 

A3 t4l| 

A 13 aA u(^lP B + : 
~*!3^** • 

hi G . 61 
0 1 2 

fl £^AAR 


-4- 4 0.04644864 -4- 4 0.01327104 -A-) 4 
t 2 t. 3 - tty 

: 015120 -A. + O.ooi8l44o -4 
4 43 

96tM - 1.19047278) 


.015150 4 0.17372192 A 


A - 0.26173519 4 

- - 0.23433120 - 

4 H 


4 E 4 1.78466080 0.33461552 ^ 4 

■ 4 0.04409856 4 




a 13 2j 

a 13 2a 33 


2.27957512 4 0.34764800 | - 0.16840704 pj 4 
4 2.90765532 4 0.85941696 | - 0.26121856 p - 0.13808882 
p - ■ 1536tP - 1.35762816 - 0.38900736 | 4 0.13229568 pj 
^165888 G + g^.44tQ 4 0.37195776 




A 33 2 A a 


4 0.28035072 4 4 0.11260384 -A- 4 0.05142528 -A-l 

4 t 2 t 3 / 

.33634816 4 0.00372096 A 4 0.11280384 -A. 4 0.05142528 -A-^ 4 

- \ t- 4 t 2 t 3 / 

a 33 ^ff^ S 6 h 4 2.8l8o6848 4 0.49208832 A + 0.25878528 p 4 0.13834592 p 4 0.03981312 h 

.87557832 - 0.30465792 £ - 0.36683712 pj 4 

Au-Hj" ' 


4-0. 

■ 2.2516876B 4 0.42723072 A - 0.0993168 o A + o. 

\ •*- 4 t 2 

a 33 2 a 13^^5S o 4 3.11637312 4 


,(*§“■ 


.42723072 A - 0.0993168o -A + 0.05142528 4 

4 t 2 t 3 / 

.03692240 A - 0.16510608 — - 0. 08204112 4 

A 33 2 a 3 i (^ 4 ^ h + 2 -9l644928 4 0.97228800 A + 0.06590592 p 4 0.13934592 p 4 0.03981312 + 

a 11 a 13 u t^l(" 576D + 2.57143464 - 0.94285944 lj 4 A^A^o + 7 68tN - 2.32002560 - 0.34228224 ij 4 

a 11 a 33 u rfr(-^ E + 2.09204160 - 1.19231904 |j + a 13 A 31 o e + 2.28251840 - 1.15422496 

*13*33° + 2 • 360536 - 0-29949696 a - 0.49552128 4 

*31*33° F + 1#Wfia72 . i.itg^gg 1 . 0.26459236 4 

4^13*31^^^ 2 + 3.911072 - 0.40761344 A + 0.08819712 -Aj 4 
*11*13*33^1^ ? + 4.10292624 4 0.53652480 A - 0.35664768 -A.\ + 


\ t 

A u A i3 A 3i( S ^ f + 3.911072 - 0.40761344 A 

A ll A 13 A 33 ? + 4.10292624 4 

A nA3lA 3 3^Ap 0 4 4. 

A ^ A Jm & 


1.31941632 


n rrrakvaio 


rt|H 
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E + 1.42932978 + 0.55309932 | - 0.13175400 i - o. 1451448 o - 0.02992983 + 

A 3! 3 ( 66 ^ 5 ~ H + °- 6225568 ° + °- o8la28 1 z £ + 0.08626176 + 0.04644864 + 0.01327104 

A 33 3 [^^ h + !• 33485408 + 0.64514880 i - o. 04784184 i - 0.09217152 -rj - 0.01550664 
All 20 t^T^ -2680 + 1- £3571732 - 0.47142972 + Au 2 a 13 ^ 3 ^ e + 1.89610200 + 0.13562784 | - 0.2 

A H 2A 31 ^pr ' e + 1.80954432 - 0.23556096 | + 0.04409856 + 


a 11 2a 31^^^^ f + 1.80954432 - 0.23556096 | + 0.04409856 + 

a U 2a 33 F + 1-91459592 + 0.22063104 i - 0.18427592 ij + 

A i3 2,J ' : E^l(" s8aD + 2 - l8 7502l7 + 0.48905445 | - 0.39710994 ^ - 0.15087705 4jj + 
a 13 2a h(^^ e + 2.72042070 + 0.69054222 i - 0.34588152 - 0.14399928 + 

A 13 2A 3i(^j^ F + a -90765532 + 0.85941696* | - 0.26122B56 i - 0.13808882 + 

A 1 3 2 A 33 F + 4.43507589 + 1.84144320 | - 0.28275813 ^ - 0.41674986 ^ - 0.08774370 pj) + 
A 31 2 U E + 0.74723136 - 0.74563584 i - 0.13229618 + 

A 31 2 An ^ .1 6 5^88 G + 2.15970816 - 0.03670656 i + 0.11260384 -i + 0.05142528 ij + 
A 31 2 Ai 3 ^^8 0 + 2.47248768 + 0.50808576 | + O.08864l6o i + 0.05142528 + 


A 31 2 An^ 1 8 5888 c + 2. 1 597o3l6 - 0.03670656 i + 0.11260384 + 0.05142528 A-j + 

a 3 1 2 Ai 3 ^™ o + 2.47248768 + 0.50808576 | + O. 08864l6o i + 0.05142528 + 

a 31 2a 33 ^2^656 H _ 2 + 0.97228800 1 + 0.06590592 + 0.13934592 ^ + 0.03981312 + 

A 33 2 a t^l('^^ F + 1 -53456444 + 0.16303680 | - 0.48873456 ^ - if 14364378 + 

^All ^g 88 0 + 2.88453312 + 0.93,080016 i - 0.18912528 i - 0.0843048o + 

A 33 gA 13 0 + 4.58072496 + 2.01231216 | - 0.16849134 \ - 0.35067168 i - 0.08593452 + 

A 33 2a 3i(^T^ h + 3.54741120 + 1.50133824 1 + 0.00316224 i + 0.01016064 i + 0.03981312 + 

a 11 a 13 u ^^j^-576d + 2.87500632 - 0.01465512 | - 0.51822222 + 

A 11 A 31 U t^l( - ^^ E + 2 -°9204l6o - 1.19231904 + 

a 11 a 33 u t ^ tt ( - E + 2.380536 - 0.29949698 i - 0.49552128 -i\ + 


A U A 33 u H iT E + 2 -3 s °536 - 0.29949696 i - 0.49552128 + 

Aj^A^O E + 2.380536 - 0.29949696 £ - 0.49552128 + 

Aj^A^U ^5^- E + 3-77748180 + 0.65672488 | - 0.75665232 ^ - 
_ . / n . . \ 


\ - / 

Aj^A^U t^l{‘ E + 3-77748180 + 0.65672488 | - 0.75665232 ^ - 0.29417024 + 

A31A33E F + 1-75U5664 - 0.60931584 i - 0.73367424 + 

AuAijAji^^^- F + 4.10292624 + 0.53652480 - - 0.35664768 + 

A H A 13A33^^E + 5-43208248 + 1.5521241 6 | - 0.55813536 i - 0.27949212 + 

A- LL A 31 A 33 ^22i2^ o + 4.50337536 + 0.85446144 | - 0.19863360 + 0.10285056 + 

A i3A3iA 3 3 ^^p' 6 G + 5-88498624 + l.9l466l44 ^ - 0.3*^24096 i - 0.16747776 Aj 


A11A31A33I 
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TABLE 1 

FLOW PAST ELLIPTIC CYLINDER WITH b/a = 0-5 AND = 0.5 


) 

Eccentric 
angle , 

(deg) 

Velocity distribution, q/u 

Incompressible 

flow 

Variational 

method 

Hooker ’ s 
method 

Kaplan’s 
results 
(reference 2) 

0 

0 




15 

.7085 

O.6655 

O.667 

0.6683 

30 

1.1339 

1.1054 

' 1.130 

1.1310 

45 

1/3V17 

1.3731 

1.380 

1.3869 

60 

1.41H2 

1.5450 

1.513 

1.5184 

75 

1.4867 

1.6468 

1.575 

1.5809 

90 

1.5 

1.6810 

1.607 

1.5994 
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TABLE 2 

FLOW PAST AN ELLHTIC CYLINDER WITH b/a = 0.1 AND Mq = 0.8 


Eccentric 


Velocity distribution, q/U 


angle, 

n 





Incompressible 

Variational 

'Perl’s solution 

Kaplan’s 
results 
(reference 2) 

(deg) 

flow 

method 

(reference 4) 

0 

0 

r"' 

0 

0 

0 

5 

• 7244 

.6982 


.6898 

10 

.9569 

.9256 

“ 

• 9137 

15 

1.0307 

1.0026 



-.9886 

20 

1.0608 

1.0399 


1.0237 

30 

1.0839 

1.0845 

1.145 

1.0629 

4o 

1.0924 

1.1199 

1.1625 

1.0918 

50 

l;0962 

1.1526 

1.1754 

, 1.U79 

60 

1-0983 

1.1819 

1.1827 

1.1408 

TO 

1.0994 

1.2052 

1.1883 

1.1589 

80 

t 1.0999 

1.2203 

1.1903 

1.1706 

90 

1.10 

1.2255 

1.1913 

r . 

■ 1.1747 
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TABLE 3 

FLOW PAST A CIRCULAR CYLINDER WHEN Mq = 0.4 


■ 

Eccentric 


Velocity distribution, q/u 


angle, 

n 

Inc ompr e s s ible 

Present solution 

Wang's results 

Rayle igh- Janzen 
method (third 
approximat ion) 

(deg) 

flow 

with four terms 

(reference 9) 

■ 

0 

.3473 

0.2966 

0.3194 

0 

• 3104 

wlM 

.6840 

.6002 

.6349 

.6439 

■Hr. 

1.0000 

.9132 

.9410 

.9588 

4o 

1.2856 

1.2305 

1.2470 

I.2658 

50 

1.5321 

1.5381 

1.5522 

1.5613 

60 

1.7321 

1.8149 

1.8450 

1.8355 

70 

1.8794 

2.0366 

2.0974 

2.0752 

80 

1.9696 

2.1804 

2.2712 

2.2269 

90 

2.0000 

2.2303 

2.3335 

2.2840 
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Figure 1.- Elliptic coordinates. 
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Figure 2.- Velocity distribution on surface of an elliptic cylinder of 

thickness ratio 0.5 when Mq = 0.5* 
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Figure 3»- Pressure distribution on surface of elliptic cylinder of 
thickness ratio 0.5 with Mq = 0.5. 
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Figure 4.- Velocity distribution, on surface of an elliptic cylinder of 

thickness ratio 0.1 when Mq = 0.8. 
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Figure 5*- Pressure distribution on surface of elliptic cylinder of 
thickness ratio 0.1 with Mq = 0.8. 
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20 UO ^ deg 60 80 100 


Figure 6.- Flow past a circular cylinder obtained by making 
t approach when Mq = 0.4. 





